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New measurements with good statistics will make it possible to observe the time variation of cosmic 
antiprotons at 1 AU through the approaching peak of solar activity. We report a new computation 
of the interstellar antiproton spectrum expected from collisions between cosmic protons and the 
interstellar gas. This spectrum is then used as input to a steady-state drift model of solar modulation, 
in order to provide predictions for the antiproton spectrum as well as the antiproton/proton ratio at 
1 AU. Our model predicts a surprisingly large, rapid increase in the antiproton/proton ratio through 
the next solar maximum, followed by a large excursion in the ratio during the following decade. 

PACS numbers: 98.70.Sa, 96.40.Kk, 95.35.+d 
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I. INTRODUCTION 

Two factors make cosmic antiprotons of special inter- 
est now. First are new experimental results |l]-Q, espe- 
cially the abundant data from the 1997 flight of BESS 
||, and the prospect of even more data from AMS ||] 
and future flights of BESS. Second is the opportunity to 
observe with good statistics the time variation of cosmic 
antiprotons at 1 AU through a solar activity maximum. 
These observations promise to provide new insights into 
solar modulation, the process by which the expanding 
solar wind modifies the energy spectrum of cosmic rays 
that enter the heliosphere. At the same time, any excess 
of antiprotons that cannot be explained as a modulation 
effect may point towards a primary source of antiprotons 
from exotic processes @-|9| . 

Cosmic ray antiprotons are good probes of solar mod- 
ulation for several reasons. First, the dominant process 
producing antiprotons is collisions of high-energy cosmic 
protons with the interstellar gas. The energy spectrum 
of antiprotons thus produced can be computed with rea- 
sonable confidence, and we therefore have a good a priori 
knowledge of the input spectrum for solar modulation 
pO|-p^,p| . Second, the antiproton production spectrum 
has a distinct peak around 2 GeV kinetic energy because 
of the high energy threshold for antiproton production in 
collisions. This is in sharp contrast to the feature- 
less monotonic spectrum of interstellar cosmic ray pro- 
tons. Third, because protons and antiprotons differ only 
in charge sign, they are ideal for studying solar modula- 
tion effects that depend explicitly upon particle charge 
sign. 

Existing evidence for charge sign dependent modula- 
tion appears in Figure 1, which displays the ratio of 
cosmic electrons to cosmic helium observed over a 25 
year period jl(J, together with recent observations of 
the electron to proton ratio made aboard Ulysses Jl7| . 



The largest variations are associated with reversals of the 
Sun's magnetic polarity (shaded bars), which occur near 
the peak of solar activity. In 1970 and again in 1990, the 
charge ratios decreased rapidly. In 1980 the ratio jumped 
upwards. If the pattern continues, another large, rapid 
increase in the negative/positive charge ratio will occur 
through the polarity reversal expected in 2000 or 2001. 
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FIG. 1. Ratio of (top) cosmic electrons to cosmic helium 
at 1.3 GV rigidity and (bottom) cosmic electrons to cosmic 
protons at 2.5 GV rigidity. Shaded areas delimit time periods 
when the Sun's poloidal magnetic field was reversing. Positive 
and negative solar polarity refer to epochs when the magnetic 
field emerging from the Sun's north pole point respectively 
outward and inward. Data are from |l6| , |l7| |. 

In the remainder of this Letter, we present a new cal- 
culation of the interstellar antiproton spectrum. This 
spectrum is then used as input to a drift model of solar 
modulation. Finally, we present our prediction for the 
evolution of the antiproton spectrum and the antipro- 
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ton/proton ratio through the upcoming solar maximum. 



II. INTERSTELLAR ANTIPROTON SPECTRUM 

In the framework of the standard leaky box model the 
continuity equation describing secondary antiproton pro- 
duction |Q can be written as 



f c 

-Jp(Ep) + —Jp(Ep) = 

. Xi 47r(m) 



Q(E P ;J P (E P )), (1) 



where A e is the characteristic escape length, Jp(Ep) de- 
notes the antiproton intensity, and is the interaction 
length for inelastic collisions of antiprotons with the in- 
terstellar gas (annihilation plus non-annihilation). The 
mean free path length is Xi{E p ) = (to) / (<7™ el (E p )) , where 

(to) and (ap 11 (Ep)) denote the target mass and inelastic 
cross section averaged over the composition of the inter- 
stellar gas, respectively. The mean escape length A e is 
taken from the recent fit to ratios of secondary to pri- 
mary nuclei by Webber et al. |l8[] . 

The source term Q is split into two parts jl3| ] 
Q{Ep\Jp(Ep)) — Q P rod{Ep) + Q S catt(Ep). Here, Q pro d 
is the source function for the production of antiprotons 
due to collisions of primary cosmic rays with the inter- 
stellar gas 
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and Qscatt takes the inelastic scattering of antiprotons on 
the interstellar gas into account 



W scatt (Ep) 
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Jp{E)dE . (3) 



The index i sums over primary cosmic ray particles (pro- 
tons and alpha-particles in our calculation) and j runs 
over all interstellar gas target particle species (H, He, C, 
N, and O). The particle abundances rij with ■ rij = 1 
are taken from the data compiled in [ p"9[ . The antiproton 
production and inelastic scattering cross sections have 
been calculated with a new version of the Dtunuc Monte 
Carlo event generator |2(],|2l| which uses Phojet Q for 
the simulation of elementary nucleon-nucleon collisions. 



III. MODULATION CODE 

The effect of gradient and curvature drifts on solar 
modulation has been intensively studied over the past 
25 years p3|-pq|. Drifts in principle can provide a nat- 
ural explanation for charge sign dependent modulation 
effects 27 1, because particles with opposite charge drift 
in opposite directions. However, in recent years there has 



been an emerging consensus that drifts may be important 
for modulation during low solar activity, but that they 
become unimportant for several years around solar max- 
imum, owing to the disordered magnetic structure of the 
heliosphere at that time p8| , p9[ . 

Recent work |3(J has challenged the conventional wis- 
dom that drifts can be ignored during high solar activity. 
This work finds that drifts produce a strong differentia- 
tion between modulation of positive and negative charges 
even during high solar activity. There may be a brief in- 
terval during the polarity reversal when the heliosphere 
is in a "no drift" state, but the approach to and through 
this state is abrupt. The observational evidence displayed 
in Figure 1 favors this point of view decisively. Indeed, 
the largest variation of charge ratios occurs during peak 
solar activity in association with the polarity reversal. 

The principal factors governing solar modulation are 
solar wind speed, tilt of the heliospheric current sheet, 
and the cosmic ray diffusion tensor (which also embodies 
the drift effect in its off-diagonal terms). For wind speed, 
we use a simple latitude dependent model consistent with 
the average properties of the solar wind |3(]] . 

The heliospheric current sheet is the surface in the so- 
lar wind that separates opposing magnetic polarities. It 
is essentially always tilted with respect to the solar equa- 
tor, by an angle that varies from about 10° at sunspot 
minimum to more than 70° during high solar activity. 
The combination of solar rotation and radial expansion 
distorts the current sheet into a wavy "ballerina skirt" 
pattern, and this complex field pattern is one of the fac- 
tors that leads to greater modulation during high solar 
activity pj]| . 

The diffusion tensor K for a coordinate system with 
one axis parallel to the background magnetic field, B D , 
and the other two perpendicular to it, is 



«H 

Kj_ !polar KT 
-Ky K_L,r 



(4) 



with K\\ the diffusion coefficient describing the diffusion 
parallel to B G , Kj_, po iar and k±^4> the diffusion coefficients 
describing the diffusion perpendicular to the background 
magnetic field in the polar and radial/azimuthal direc- 
tions respectively, and kt the drift coefficient. We use 
quasi-linear theory and a slab/two-dimensional geometry 
for the turbulence to calculate «|| |32| and recent theoreti- 
cal results |3^] for the spatial variation of the fluctuations 
in the field as well as its correlation length, both quan- 
tities which appear in Kit. The perpendicular diffusion 
coefficients, K^ iPO i ar and Kx,r<f>t are assumed to have the 
same spatial dependence as K||, but a different rigidity 
(R) dependence. For the drift coefficient kt we assume 
that kt I particle speed is proportional to R 3 at low rigid- 
ity, rolling over to R 1 at high rigidity. A detailed discus- 
sion of a diffusion tensor similar to the one used in this 
work can be found elsewhere ]3C| ]. 

Diffusion parameters were determined by fitting model 
results to proton observations at 1 AU. This procedure 
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starts from an assumed interstellar proton spectrum. Be- 
cause a series of recent results |34| |38[1 indicates that the 
proton spectrum in the range 10 to 100 GeV/nucleon 
is significantly lower than previously assumed |}9], we 
have correspondingly revised downward the standard in- 
terstellar proton spectrum | flo[ in this energy range. Our 
assumed interstellar hydrogen spectrum, shown by the 
solid line in Fig. 2a, fits smoothly to the original result 
of jh)| in the low energy region and to the data 34 3^] 
at high energy. 

The fit to the solar minimum proton data p2| shown 
in Figure 2a, is obtained with a tilt angle of 10°. To 
model spectra near solar maximum, a tilt angle of 70° 
is used [4^] while «||, n± tPO i al and 

fc-L,r(p £1X6 all reduced 
by a factor of 2/3. Only the sign of the magnetic field is 
changed when going from one solar polarity epoch to the 
other, i.e. it is only the direction in which particles drift 
that distinguishes the two polarity states in the current 
study. 



IV. RESULTS 

In Fig. 2b the prediction for the local interstellar 
antiproton intensity is shown. Our calculation is in 
good agreement with the results reported in [fL3|| , with 
a slightly higher value below one GeV, and it agrees very 
well with the calculation of || at all energies. We show 
here only the central value of our calculation. The theo- 
retical uncertainties due to the limited experimental data 
available for several quantities entering the antiproton 
intensity calculation will be discussed elsewhere [}43fl , as 
well the small differences with other calculations. The 
antiproton spectra at 1AU are also shown in Fig. 2b 
with the same coding as in Fig. 2a, and Fig. 2c shows 
the corresponding p/p ratios. 

Figure 3 displays the predicted dependence of the pro- 
ton and antiproton intensity at 1 AU (relative to inter- 
stellar level) upon tilt angle of the heliospheric current 
sheet, as well as the predicted dependence of their ra- 
tio. The abscissa values have been arranged so that the 
curves have the appearance of two successive solar cycles 
evolving in time. At the solar maximum of 2000, the solar 
polarity switches from positive to negative. The upper 
panel displays a well known feature of drift models [£4| : 
the curves for positive charges are broad during epochs of 
positive solar polarity (1990's), and pointy during epochs 
of negative polarity (decade beginning in 2000). The op- 
posite relationship holds for negative charges. 

Another difference is that protons have a greater mod- 
ulation amplitude (~ 4x between solar maximum and 
solar minimum) than do antiprotons (~ 2x) E5J. This 
stems from the differing character of their input spec- 
tra: The interstellar antiproton spectrum, with its peak 
around 2 GeV, is "hard" compared to the proton spec- 
trum which has many particles below 2 GeV. This feature 
p5[ has been suggested by as a way to distinguish a 



component of primary antiprotons from an exotic source 
(such as annihilation of neutralinos in the galactic halo or 
evaporation from primordial black holes) from the galac- 
tic secondary component discussed in this paper. Real- 
ization of this idea depends on the extent to which the hy- 
pothetical primary spectra are softer than the secondary 
spectra. Recent calculations of the spectra from neu- 
tralino annihilation [^],^| show progressively harder spec- 
tra, making this signature correspondingly less helpful. 
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FIG. 3. Predicted dependence of (top) the proton and an- 
tiproton intensity and (bottom) the antiproton/proton ratio 
at 1 AU upon the tilt angle of the heliospheric current sheet 
for 1 GeV kinetic energy. Intensities are relative to inter- 
stellar level. Abscissa values are arranged so that the curves 
mimic the expected time variation through two solar cycles 
of opposite magnetic polarity (1990 to 2010). 



V. SUMMARY 

Secondary galactic antiprotons provide a powerful 
probe of solar modulation. Protons and antiprotons have 
sharply different interstellar spectra. They also drift in 
opposite directions because of their opposite charge signs. 
The combination of these effects implies that the antipro- 
ton/proton ratio should display a much more interest- 
ing evolution during the next 10 years than it did dur- 
ing the 1990's, when the ratio was nearly constant. As 
we proceed through the sunspot maximum and polarity 
switch expected about 2000, we predict that this ratio 
will rapidly increase by a factor of about 3. Then, dur- 
ing the following decade, it will display a large excursion 
closely tied to the variation of the current sheet tilt an- 
gle. Actual observation of these variations would be a 
stunning validation of the importance of drift effects in 
solar modulation at all phases of the solar activity cycle. 
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